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Abstract. In this paper we study the combined effects of size and shape of small solid state particles on the absorption, 
emission and scattering characteristics. We use the statistical approach to calculate these optical properties. In this approach the 
average optical properties of an ensemble of particles in random orientation are represented by the average optical properties 
of an ensemble of simple shapes. The validity of this approach is studied in detail for a uniform distribution of hollow spheres 
where the fractional volume of the central inclusion is varied. We apply the results to two different areas of interest, i) infrared 
spectroscopy and ii) polarization of scattered light. The effects of particle size and shape on the optical characteristics are 
discussed. We compare the results using the distribution of hollow spheres with those obtained by using randomly oriented 
spheroids. Also we compare the results with observations and laboratory measurements. The distribution of hollow spheres is 
very successful in reproducing laboratory measurements of the scattering angle distribution of the degree of linear polarization 
for incident unpolarized light of randomly oriented irregular quartz particles. Furthermore, we show that we are able to derive 
the size distribution of dust grains by fitting the measured degree of linear polarization using computational result for hollow 
spheres. It is shown that the distribution of hollow spheres is a powerful tool for studying light scattering, absorption and 
emission by cosmic dust grains and in particular when large numbers of particle parameters need to be considered since the 
computational demand of the distribution of hollow spheres is small. 


1. Introduction 


Dust is an important constituent of many astrophysical as well 
as atmospheric environments. The interpretation of observed 
infrared spectra, images and polarization requires knowledge 
of the scattering, absorption and emission properties of dust 
particles of various compositions, shapes and sizes. Optical 
properties calculated with the widely used Mie theory, which 
is valid for homogeneous spherical particles, often do not give 
satisfactory agreement with measurements and observations. 
Therefore, complicated methods have been developed to cal¬ 
culate the optical properties of irregularly shaped particles such 
as the Discrete D ipole Approximation (PDA) and the T-matrix 
method (see e.g. Purcell & Pennvpackeil 197^ Drain^ll988t 


iMishchenko et al.i boool l2002lTwriedl l20Q2t iKahnertl I 2 OO 3 I) . 

For realistically shaped particles these methods are all very 
computationally expensive. Another approach to obtain opti¬ 
cal properties of ensembles of irregularly shaped particles is 
the so-called statistical approach. In this approach the average 
optical properties of an ensemble of irregularly shaped parti¬ 
cles are represented by the average properties of a distribu¬ 
tion of simply shaped particles with a given shape distribu¬ 
tion. Examples of simple shapes used in this way are spheroids 
and circular and hexagonal cylinders dMishchenko & Travis 
1994 IMishchenko et all 1996at[Kahnert et al.ll2002tlKahnerl 


2004]) . In iMin et alJ ( l2003ir it was shown that the statistical 


approach is very successful in calculating the absorption cross 
sections of an ensemble of randomly oriented very small parti¬ 
cles, i.e. in the Rayleigh regime, as a function of wavelength. In 
that paper also results of calculations were presented for a uni¬ 
form distribution of hollow spherical particles, varying the vol¬ 
ume fraction occupied by the central inclusion. This shape dis¬ 
tribution reproduces the measured absorption cross sections of 
small crystalline forsterite particles remarkably well. Because 
of their spherical symmetry the optical properties of hollow 
spheres can be obtained easily by using an extension of Mie 
theory. This makes this grain model an ideal candidate for stud¬ 
ies where a large number of calculations for various grain sizes 
and/or compositions is required. 

Ijonesllfl^ used hollow spherical grains to model grain 
porosity in the ultraviolet to near infrared. The computational 
results using hollow spheres were compared with those ob¬ 
tained by using effective medium theories. The main motiva¬ 
tion of that study was to find out if porous grains could enhance 
the interstellar extinction per unit mass. This turned out to be of 
much astrophysical significance in view of the s o-called ’car¬ 
bon cr isis’ which was brought up 7 years later bv ISnow & Witll 
dl 99.51) . It was found that, although by using porous particles 
the extinction is enhanced at some wavelengths, it is decreased 
at other wavelengths, consist ent with the expectation from 
the Kr amers-Kronig relations llPurcelll[l969l) . However. Lion e si 
(H) argues that this problem might be solved by changing 
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the size distribution of interstellar grains. We will discuss the 
effect of grain shape and porosity on the interstellar silicate 
extinction feature in the mid infrared by applying the hollow 
sphere model. 

In this paper we will use a distribution of hollow spheres 
to calculate absorption, extinction and scattering properties for 
various particle sizes beyond the Rayleigh regime. We com¬ 
pare our results to those obtained by using a shape distribution 
of spheroidal particles and to laboratory measurements of the 
scattering properties of irregular particles. 

In Sect. 12 we introduce the dust shape models employed 
in this paper. Then we will concentrate on two main areas; i) 
absorption and extinction spectra (Sect. |3l and ii) the degree 
of linear polarization (Sect.^J. Some examples of applications 
in both areas are provided. In Sect. |3 we discuss the physical 
interpretation of the results. 


2. Dust shape models 

In order to understand the effects of particle shape and inter¬ 
nal structure on the emissivities of dust grains we consider two 
simple dust models. As a first approximation we change the 
internal structure of homogeneous spherical dust grains in the 
most simple way by using a distribution of hollow spherical 
particles, varying the volume fraction occupied by the central 
inclusion. In the second distribution we change the external 
shape of the homogeneous spherical particles by using a dis¬ 
tribution of spheroidal particles, varying the aspect ratio. In 
these shape distributions the material volume of the particles is 
kept constant. Both shape distributions have been studied in the 
Rayleigh limit, i.e. when the particles are very small compared 
to the w avelength inside and outside the particle. In iMin et alJ 
<2003ai) it was found that in this limit the absorption properties 
for the distribution of spheroidal particles and for the distri¬ 
bution of hollow spherical particles are very similar and that 
both shape distributions provide good agreement with labora¬ 
tory measurements of the absorption cross sections as a func¬ 
tion of wavelength in the infrared part of the spectrum. This 
suggests that the exact shape of the particles is to a certain de¬ 
gree not important as long as we break the perfect symmetry of 
a homogeneous sphere. When applying the statistical approach 
one needs to realize that the particle shapes used contain only 
very limited information on the real shapes of the dust grains 
and that additional research has to be done in order to decode 
this informat ion into physical cha racteristics of the dust grains 
(see also e.g. iKahnert et al.Ll2002l) . 


2.1. Distribution of Hoiiow Spheres (DHS) 


The simplest internal deviation from a homogeneous sphere 
is a hollow spher ical shell with the same material volume. In 
Min et al. 1 200.3 ah a Distribution of Hollow Spheres (DHS) is 


used to calculate the absorption cross sections for very small 
particles. In this shape distribution we simply average over the 
fraction / occupied by the central vacuum inclusion in the par¬ 
ticle from zero to some value /max ^ / giving equal weights to 
all values of / while the material volume of the particle is kept 


the same. This means that the particles with higher values of / 
have a larger outer radius. Thus, the distribution is given by 

( l//max, 0 < / < /max, 

«(/) = j (1) 

(0, / > /max, 

where n(f)df is the fraction of the number of particles in the 
distribution between / and / -H df. Unless stated otherwise, in 
this paper we will take /max = 1, the most extreme shape dis¬ 
tribution. We have to note here that particles with f - I and 
a finite volume have an infinitely large outer radius. For parti¬ 
cles much smaller than the wavelength, integrating the optical 
proper ties up to / = 1 can be done analytically jMin et al.L 
l2003al) . For larger particles numerical computations are neces¬ 
sary, but these are not possible for / = 1. However, we then 
chose /max large enough to reach convergence to the values for 
/max = 1 ■ For most cases we considered, it suffices to integrate 
up to / = 0.98. The choice of the shape distribution we take is 
motivated by the fact that for many applications one of the pa¬ 
rameters of interest is the dust mass. Therefore, it is convenient 
to choose a shape distribution in which the total mass of the 
particles is conserved. Other choices of the shape distribution, 
e.g. conserving the surface area of the particles, might be con¬ 
venient in other cases. The equations for the optical cross sec ¬ 
ti ons in the Rayleigh limit are provided bv lvan de Hulstlill9.37n . 
In lAden & Kerkeif(ll95ll) the general equations for light scat¬ 
tering and absorption by layered spheres of a rbitrary size are 
derived as a simple extension of Mie theory (iMiei 119081) . By 
taking the refractive index of the core of a layered particle to 
be that of vacuum these general equations can be used to cal¬ 
culate the optical properties of hollow spherical particles for 
almost all particle sizes. 

The calculations were done using the l ayered sphere code 
for wh ich the basic ideas are explained in Ploon & AckermanI 
EUL We successfully perf ormed the Cloude test (see 
iHpyenier & yan der Med I2OOOI) on the 4x4 scattering ma¬ 
trix obtained with this code. Also we checked that the com¬ 
puted yalues of the scattering matrix obtained for forward and 
backward scattering are in agreement with theoretical yalues 
deriyed from symmetry arguments. In addition, we considered 
the limiting case / —> 0 and compared this with calculations 
for homogeneous spheres, all with satisfying results. 

2.2. Uniform Distribution of Spheroids (UDS) 

The simplest external deyiation from a homogeneous sphere 
is a homogeneous s pheroid with the same material yolume. 
In iMin et al.l il2003bl) a method was deyeloped to calculate the 
optical cross sections of a Uniform Distribution of Spheroids 
(UDS) for eyer y particle size using a com bination of the T- 
matrix method iMishchenko et all Il996bh and some approx¬ 
imate methods. In the UDS we ayerage oyer the aspect ra¬ 
tios (i.e. the ratio between the major and minor axis) of 
the spheroids giying equal weights to all oblate and prolate 
spheroidal shapes. The material yolume of all parti cles in the 
shape dis tribution is the same. For details we refer to lMin et"^ 
il2003allJ) . Although the ayerage absorption cross sections for 
this distribution can be calculated, the computations become 
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x = 0.1 x=1.0 x=2.0 x=10.0 



Re(m) Re(m) Re(m) Re(m) 


Fig. 1. The shape averaged absorption efficiency as a function of the real and imaginary part of the refractive index for vari¬ 
ous values of the dimensionless size parameter x - 2nryjX, where ry is the radius of a volume equivalent sphere and A is the 
wavelength of incident radiation. The upper panels are for homogeneous spheres, the lower panels for the distribution of hollow 
spheres (DHS). 


very slow for intermediate sized particles (with sizes compara¬ 
ble to the wavelength of radiation) where we still have to apply 
the T-matrix method. 

In this paper we will average the spheroidal shapes over all 
prolate and oblate shapes up to an aspect ratio of 10. This is 
at the moment the most extreme shape distribution for which 
calculations can be done in a reasonable time. For the T- 
matrix calculations we used the well tested code developed by 
iMishchenko & Travis! ( Il998l) . 


3. The absorption, extinction and emission spectra 

Absorption, extinction and emission spectra of dust provide 
crucial information on the size and composition of the dust 
grains. The thermal emission spectrum of a dust grain in ther¬ 
mal equilibrium averaged over all orientations is proportional 
to the orientation averaged emission cross section times a 
Planck function. From Kirchhoff’s law we know that this emis¬ 
sion cross section is equal to the orientation averaged absorp¬ 
tion cross section. The size of a particle, ry, is defined to be the 
radius of a sphere with the same material volume as the parti¬ 
cle. We define the average absorption cross section, (Cabs), to 
be the shape- and orientation-averaged absorption cross section 
of a collection of particles with a shape distribution and a fixed 
value of the material volume equivalent radius. In this paper we 
will always first average over all particle orientations. After this 
we will average over particle shape, which is denoted by one 


bracket (..). For calculations of the optical characteristics of a 
particle it is useful to define the dimensionless size parameter 



( 2 ) 


where A is the wavelength of incident radiation. The mass ab¬ 
sorption coefficient k is the average absorption cross section per 
unit mass 


K — 


(Cabs) 

M 


(3) 


Here M is the mass of the particles in the collection with the 
same ry. 

The orientation averaged absorption efficiency of a parti¬ 
cle defined as Qabs = Cabs/G^ with G the orientation averaged 
geometrical shadow, only depends on the shape, size parame¬ 
ter and complex refractive index m of the particle. Because the 
distributi on of normals to the surface of all convex particles is 
uniform dvan de Hulsllll95% . the absorption efficiency of very 
large convex solid particles averaged over all orientations is the 
same as for spheres and can be written as dBohren & Huffm^ 
Il983h 


eabs = l-;^, (4) 

where "R is a reflection factor which does not depend on the 
size and shape of the particle but only on the complex refrac¬ 
tive index m. Eq. 0 also holds for particles with one or more 
inclusions when the outer surface of the particle is convex (e.g. 
hollow spheres). Therefore, for very large particles the absorp¬ 
tion efficiency and cross section averaged over all orientations 
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Fig. 2. Absorption cross section per unit mass for amorphous 
olivine for different shape distributions and grain sizes. The 
solid line is for the distribution of hollow spheres (DHS), the 
dotted line for the uniform distribution of spheroids (UDS) and 
the dashed line for homogeneous spheres. 


as a function of wavelength will have the same behavior for all 
particles with a convex outer surface. 

In Fig. n the shape averaged absorption efficiency, (Qabs); 
is plotted as a function of the real and imaginary part of the 
refractive index for various values of the size parameter for ho¬ 
mogeneous spheres (upper panels) and for the distribution of 
hollow spheres (DHS, lower panels). We clearly see that there 
is a large difference between the contours for homogeneous 
spheres and those for the DHS when x is small and that this 
difference slightly disappears when we go to larger values of x, 
which is to be expected from Eq. 0 . The difference between 
the two types of contours is already small for x = 10, which is 
much too small to expect Eq. 0 to be valid. Based on the as¬ 
sumptions used to derive Eq. 0 , it is expected to be valid only 
when X > 60. This suggests that for particles with a size param¬ 
eter larger than about 10 with a convex outer surface the behav¬ 
ior of the absorption cross section as a function of wavelength 
is to a certain degree independent of particle shape, which al¬ 
lows us to use homogeneous spheres to calculate the absorption 
cross sections for x > 10. 


3.1. Example: amorphous silicates 


As a first example we consider the infrared mass absorp¬ 
tion coefficient of amorphous olivine (Mg 23 ,Ee 2 - 23 ,Si 04 with 
0 < y < 1) which is one of the most abundant species found in 
astronomical environments. This silicate is observed by look¬ 
ing at the thermal emission from, for exam ple, comets (see 
e.g. lCrovisier etj lll997HHanneretal.Hl994, proto-planetar y 


disks (see e.g. lBouwmaneta ll l200lUvan Boekel et aLLboOSl) 

asymptotic giant bran ch stars (see e.g. iMolstereTd . Il999l: 
iHoogzaad et all l2002h as well as by l ooking at the extinc- 
tion by interstellar dust grains (see e.g. IPraine & Leel 1 1984 



DHS (y=1.0) 


DHS (y=0.5) 


UDS (y=1.0) 


UDS (y=0.5) 


Homogeneous 
Spheres (y=1.0) 


Homogeneous 
Spheres (y=0.5) 


Fig. 3. The normalized, continuum-subtracted mass extinction 
coefficients for yarious grain shape distributions and chemical 
compositions (black curyes). The grain size ry = 0.1 /rm. Also 
plotted is the normalized, continuum subtracted extinction of 
interstellar dust grains as obsery ed towards the Galactic center 
taken from fKemner et al' ] (12004 (gray curyes). Th e refractiye 
indice s as functions of wayelength are taken from lager et ^ 
(I 2 OO 3 I) for the y = 1.0 data and from lPorschner et al.l(ll995h fo 
the y = 0.5 data. The continuum that was subtracted from all 
spectra is taken to be a straight line from 8 to 13 fim. The nor¬ 
malization is done such that the maximum yalue of all curyes 
equals unity. The curyes are shifted with respect to each other 
with interyals of 0.5. 


lEi & Drainell200lUKemner et al.Ll2004 . The exact behayior of 
the mass absorption and extinction coefficients as functions of 
the wayelength for yarious grain sizes is therefore crucial in un¬ 
derstanding the formation and processing of grains in these en- 
yironments. The obseryed spectrum is dominated by two broad 
features, one at 9.8 ;um and one at 17.5 yum. The 9.8 yum fea¬ 
ture has been studied yery extensiyely because the atmospheric 
window around lOyum allows ground based obseryations. 


Fig El shows the mass absorption coefficient, k, as a func¬ 
tion of wayelength for amorphous oliyine grains with yari¬ 
ous sizes and y = 0.5. Th e refractiye indices as a function of 
wayele ngth are taken from lPorschner et~an lll995l) : ljager et all 
il2nn3h . We see that for the smaller grain sizes, the 17.5yum 
feature broadens significantly when going from homogeneous 
spheres towards the distribution of hollow or spheroidal parti- 
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cles. This becomes important when identifying features in the 
region 20 - 30 jum. 

It should be noted that all spectroscopic features that we 
considered are broadened with respect to those calculated using 
homogeneous spherical particles. Actually, it might be argued 
that the correct way to put it is that the features as calculated us¬ 
ing homogeneous spheres are sharpened with respect to those 
of all other particle shapes. The reason for this is that small ho¬ 
mogeneous spherical particles exhibit a rather artihcial sharp 
resonance close to when the refractive index m - V^. For 
irregular particles the mass absorption coefficient is a smoother 
function of the refractive index which has its maximum approx¬ 
imately when the bulk refractive index is at its maximum. This 
causes the difference in feature shape observed in, for example, 
the 18 micron feature in Fig.|2 

Many sources show a 9.8 fim extinction or emission fea¬ 
ture that is remarkably compatible with amorphous olivine 
grains that are very small (ry < 0.1//m), homogeneous in 
composition, spherical in shape and have a magnesium/iron 
content with y = 0.5. When we change the shape of the par¬ 
ticles, the 9.8 jum feature shifts towards longer wavelengths 
and broadens (see Fig. □l which is not compatible with ob¬ 
servations. This is puzzling since for example images of in¬ 
terplanetary dus t particles show tha t these grains have a very 
irregular shape il Warren et al .Hi 9941) . In addition grains in in¬ 
terstellar space produce observable linear polarization in the 
forward scattering direction which can only be explained us¬ 
ing non-sphe rical dust grains that are aligne d to a certain de¬ 
gree (see e.g. lHildebrand & Dragovanlll995li . Fig.|3shows the 
normalized continuum subtracted mass extinction coefficient in 
the 10//m region of amorphous olivine grains with y = 0.5 and 
y = 1.0. Note that here we plot the mass extinction coefficient, 
which is just the average extinction cross section per unit mass. 
For very small particles the scattering is negligible, in which 
case the mass extinction coefficient is equal to the mass absorp¬ 
tion coefficient, k. The figure shows the normalized continuum 
subtracted mass extinction coefficient for small homogeneous 
spherical particles, for the uniform distribution of spheroids 
(UDS) and for the distribution of hollow spheres (DHS) all with 
ry - 0.1 fim. The continuum that was subtracted is taken to be 
a straight line between 8 and 13 //m. The normalization of all 
curves in Fig. 0 is such that the maximum value equals unity. 
Also plotted in this figure is the normalized c ontinuum sub¬ 
tracte d extinction towards the Galactic center jKemper et al.L 
I 2 OO 4 I) . We see that the extinction spectrum towards the galac¬ 
tic center is fairly well reproduced by both curves with homo¬ 
geneous spherical particles. The feature obtained by using the 
spheroids or hollow spheres with y = 0.5 is much broader than 
the galactic extinction feature. Also the maximum for these 
grains is shifted towards longer wavelengths than observed. For 
the magnesium silicates (y = 1.0), the agreement with the in¬ 
terstellar extinction feature for the DHS and UDS shape dis¬ 
tributions is better. Because of the observed degree of linear 
polariza tion the interstellar grai ns are most likely non-spherical 
(see e.g. lGreenberg & 0.119961) . A possible explanation for the 
shape of the observed interstellar extinction feature might be 
that the grains in the interstellar medium have a very high mag¬ 
nesium over iron ratio. Note that the particle ensembles con¬ 



X(/um) X(/U.m) 


Fig. 4. Mass absorption coefficient, k, for various sizes of crys¬ 
talline forsterite particles. The solid line is for the distribution 
of hollow spheres (DHS), the dotted line is for the uniform dis¬ 
tribution of spheroids (UDS ). The refractive index as a function 
of wavelength is taken from IServoin & Pirioul(ll973h . 


sidered here are in random orientation and would therefore not 
produce any linear polarization in the forward scattering direc¬ 
tion. However, aligned, non-spherical dust grains produce an 
extinction spectrum which is more similar to that of randomly 
oriented non-spherical or non-homogeneous grains than to that 
of homogeneous spheres. 

3.2. Example: crystalline silicates 

Another example of the importance of grain shape for the ab¬ 
sorption and extinction spectra are the spectra of crystalline 
silicate particles. Small crystalline silicate grains have a very 
clear spectral signature in the infrared part of the spectrum 
which makes the m easy to detect in for e xample thermal emis¬ 
sion spectra (e.g. iBouwman et d I I 2 OO 3 I and references at the 
start of Sect. HD Crystalline silicates can form only under cer¬ 
tain conditions from gas phase condensation or by annealing of 
amorphous silicates. The abundance of crystalline silicates is 
therefore a good probe for the processing history of the dust 
grains. 

The mass absorption coefficient as a function of wave¬ 
length for small crystalline silicate grains shows many very 
narrow and strong resonances in the infrared part of the spec¬ 
trum. Because of the strength of the resonances, their shape 
and wavelength positions are very sensitive to the particle 
shape. For example, calculations of the mass absorption coef¬ 
ficient as a function of wavelength using homogeneous spher¬ 
ical particles completely fail to reproduce laboratory measure¬ 
ments on small cr ystalline forsterite grains JFabian et alll2001t 
iMin et aklEnn^ . A yery successful shape model which is fre¬ 
quently used to calculate the absorption spectra of these grains 

Distribution of Ellipsoids (CDE) 
This distribution takes all possible 


is the so-called Continuous 
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Fig. 5. Same as Fig.E]but only for /I = 8 - 13.5 yum. 


ellipsoidal shapes into account. Although the measured absorp¬ 
tion spectrum is very well reproduced using this shape distri¬ 
bution, a major disadvantage is that calculations for this distri¬ 
bution are restricted to the Rayleigh domain since it uses ellip¬ 
soidal shapes, and it is thus not suitable for studying combined 
size/shape effects. We pe rformed calcula tions using the shape 
distribution suggested bylu et al.lll2002i) . However, the results 
for, for example, the wavelength positions of the crystalline sil¬ 
icate resonances were poorly reproduced, in contrast with using 
CDE. The shape distributions we use h ere also reproduce the 
measurements on small forsterite grains ilMin et al.[l2003al) and 
can be used outside the Rayleigh domain. 

Fig. 13 shows the mass absorption coefficient as a function 
of wavelength for various grain sizes. The refractive index as a 
function of wavelength is taken from IServoin & Pirioullll973h . 
The solid curve is for the distribution of hollow spherical par¬ 
ticles (DHS), the dotted curve is for the uniform distribution 
of spheroids (UDS). We see that the strong resonances visible 
for very small particles disappear when we increase the parti¬ 
cle size. We notice that the resonances at shorter wavelengths 
disappear first. Also, the stronger resonances disappear before 
the weaker resonances do. These effects are due to the fact that 
the important parameter which determines the absorption prop¬ 
erties is \2n m ry jA\. So the strength of the resonances is influ¬ 
enced by the wavelength, the refractive index and the size of 
the particle. This effect causes the ratios between the various 
resonances to change when we increase the particle size. In the 
interpretation of observations this ratio can be used as an indi¬ 
cator of the typical particle size. 

We now compare the curves calculated for the uniform dis¬ 
tribution of spheroids with those calculated for the distribution 
of hollow spheres in Fig. |3 First we note that for very small 
particles (ry = 0.1 jum) the differences are very small. When 
the grain size is increased the trends seen in the two curves 
are the same, as is described above. However, we see that the 
resonances for a given grain size obtained from the distribu¬ 
tion of hollow spheres are generally stronger. Since in the DHS 


we average over all possible inclusion fractions from 0 to 1, 
we always include spheres with extremely thin shells. These 
thin shells will, to a certain degree, simulate the optical proper¬ 
ties of very small homogeneous particles, resulting in stronger 
spectral features. This can be understood by considering the 
particle as a combination of small volume elements, each inter¬ 
acting with the electro magnetic field like dipoles (the Discrete 
Dipole Approximationi lPurcell & PennvnackeilllQTsI) . If these 
volume elements are close together, the interaction between 
the volume elements is important. However, if the volume el¬ 
ements are distributed over a larger area, as in a thin shell, the 
interactions between the volume elements are less important 
and the total interaction will be closer to the sum of the sepa¬ 
rate dipoles. In a thin shell the volume elements are distributed 
over a larger volume area than in a homogeneous spheroid with 
the same material volume, so the interaction of the light with 
the separate dipoles becomes more important. 

As already mentioned, the 10 micron region is very impor¬ 
tant in astrophysics because of the possibility of ground based 
observations. Fig. 0 is a blowup of the 10 micron region. It 
shows that for the detailed shape of the spectrum, there are dif¬ 
ferences between the UDS and DHS results. In both sets of 
curves the trends as discussed above are still clear. Besides 
a change in the feature strengths and ratios, we also observe 
a shift in the wavelength positions of the resonances. For ex¬ 
ample the 11.3 jim feature shifts towards 11.45 - 11.5 jum (de¬ 
pending on the shape distribution used) when the grain size 
is increased from ry = 0.1 //m to ry = 1.0 pm. A similar shift 
from 11.3 to 11.44pm is observed when the chemical composi¬ 
tion is changed from magnesium rich (forsterite, v = 1) to iro n 
rich (fayalite, y = 0) jFabian et all 1200 It Honda et al. I20n4ll . 
Therefore, it is necessary to consider multiple resonances, at a 
large range of wavelengths, to disentangle the effects of par¬ 
ticle size and chemical composition. Because we are looking 
at relatively short wavelengths, the spectral diagnostics of the 
emission resonances of forsterite are already reduced signif¬ 
icantly when the volume equivalent radius of the particles is 
~ 2 pm. Therefore, the mineralogy that can be deduced from 
10 pm spectroscopy is limited to particles smaller than a few 
micro n. This has to be kep t in mind when analyzing these spec¬ 
tra (cf. lHonda et aklEoM) . 


4. The degree of linear polarization 

Another area we wish to consider is the degree of linear po¬ 
larization of the light scattered by dust grains in the visible 
part of the spectrum, which is a frequently used observable. 
For comets it is possible to obtain measurements for various 
angles of scattering by doing measurements at different orbital 
phases of the comet. For discussions on the polarization of 
light scattered by cometary du st and the possible implica tions 
for the grain properties s ee e.g. ll .evasseur- Regour(j i ll999l) and 
iKolokolova et al Jll200lh . The degree of linear polarization as a 
function of scattering angle can also be measured directly in the 
lab oratory for collections of randomly oriente d particles (see 
e.g. iHoyenier et al 1 120031 IVolten et all l200ll) . In most cases 
the shape of the linear polarization curye cannot be reproduced 
by using homogeneous spherical particles. Therefore it is yery 
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Fig. 6. The degree of linear polarization of solid and hollow spheres as a function of the scattering angle and the effective size 
parameter for different values of /max- The refractive index of the m aterial m = 1.76 + 0.0241 which is typical for amorphous 
silicate with 30% iron at A - 0.65 fim (iHenning & StognienkolfloQ^) . Upper left /max = 0 (homogeneous spheres), upper right 
/max = 0.2, lower left /max = 0.5 and lower right /max = 1.0. A gamma size distribution with Veff = 0.1 was used (see Eq. ( 0 ). 
Solid and dotted curves represent the areas with positive and negative polarization respectively. 


interesting to calculate the polarization properties of inhomo¬ 
geneous or non-spherical particles. Many studies concerned 
the linear polarization of collections of spheroids or cylindep 

( iMishchenko et al.ll2002HKahnert et al.ll2002tlKahnerttl2004h . 

Here, we will focus on calculations using the distribution of 
hollow spheres and compare our results with homogeneous 
spheres and laboratory measurements on irregular quartz parti¬ 
cles. 

Figs.0and0show the degree of linear polarization of light 
scattered by collections of particles for incident unpolarized 
light as a function of scattering angle, 6, and effective size pa¬ 
rameter. Fig. 0is for a refractive index typical for olivine with 
y = 0.5, Fig . 0is for a magnesium rich olivine with y - 1.0. In 
order to reduce the effects of resonances a narrow gamma size 
distribution given by 


n(x) - C ■ 3Veff)/VeiI 


— )• 

V'effJCeff / 


(5) 


was employed. In this equation n(x)dx is the number of parti¬ 
cles with size parameter Inry!A between x and x -H dx, Xeff is 
the effective size parameter, Veff is the effective variance, which 
we take to be 0.1, and c is a normalization constant. For val¬ 
ues of the effective size parameter smaller than approximately 
unity we see in Figs. 0 and 0 the well known Rayleigh pat¬ 
tern with a maximum degree of linear polarization of 100% at 
90°. For the homogeneous spheres and for the distribution of 
hollow spheres with /max = 0.2 large areas in parameter space 
produce negative polarization. There are even values of the ef¬ 
fective size parameter for which the linear polarization is nega¬ 
tive over the entire range of scattering angles. For the most ex¬ 
treme distribution of hollow spheres (with /max = 1) we see that 
the regions of negative polarization shift towards the very small 
and large scattering angles. This is consistent with a frequently 
observed negative pol arization branch at large scattering an¬ 
gles (6> 160 - 180°1 llVolten et al.lEooillMiifioz et al.H2000t 

iLevasseur-RegourS Il999t iKolokolo^ et all 2001 ). Also we 
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Fig. 7. As Fi g. Ibibut for a refractive index of the material m = 1.67 + O.OOOli which is typical for iron poor amorphous olivine at 
A - 0.65 um lIHenning & Stognienkotll996l) . 


note that the difference between the values obtained for olivine 
particles with y - 0.7 and olivine particles with y = 1.0 de¬ 
creases when we increase /max- This suggests that for very ir¬ 
regular particles the refractive index of the material is of lesser 
importance. This agrees well with the observed similarity of 
the degree of linear polarizati on as a function of scattering an¬ 
gle for various mineral types jVolten et al.H200lh . 

The contours in Figs. and 0 for /max = 10 show that the 
degree of linear polarization for this shape distribution is very 
high (up to 90% at intermediate scattering angles) even for very 
large particles. The high degree of linear polarization seen here 
might be due to the fact that although the outer radius of the 
particles is large, the thickness of the shell is still very small 
compared to the wavelength, which causes a high degree of 
linear polarization (see also Sect. 13.2> . 

4.1. Example: the polarization of small quartz particles 

Quartz, Si02, is one of the most abundant minerals in the 
Earth’s crust. It is believed to be for med when amorphou s sili¬ 
cates are annealed to form forsterite (iFabian et ail l2000l) . The 


refractive index of quartz in the visible part of the spectrum is 
very well known, which makes it an idea l test case for com- 
parison between measurements and theorv. fVolten et af] (l200lh 
measured the complete scattering matrix as a function of scat¬ 
tering angle in the range 5-173° of an ensemble of small 
irregular quartz particles in random orientation at two wave¬ 
lengths, 442 nm and 633 nm. The data of these measurements 
are available online together w ith the particle size distribution 
obtained from laser diffraction jVolten et al.[l2004l) . From these 
scattering matrices we can obtain the degree of linear polariza¬ 
tion for incident unpolarized light 


P(6) = - 


Fim 

Fuioy 


(6) 


where 0 is the scattering angle and F„i^ is the «, kth element o f 
the 4x4 scattering matrix (for details see lVolten et al.[E004t) . 
The laboratory data provide the element F 12 . However, the 
symmetry of the system yields F 12 = F 2 i. For comparison we 
also compute the degree of linear polarization using the distri¬ 
bution of hollow spheres together with the measured particle 
size distribution. For consistency with calculations elsewhere 
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Fig. 8. The degree of linear polarization of small quartz particles versus scattering angle, Q. Shown are the measurements (points) 
and calculations (curves) using the measured size distribution and the distribution of hollow spheres with /max = 0.3 (dotted 
curve), /max = 0.4 (solid curve) and /max = 0.5 (dashed curve). The left figure is for a wavelength A = 422 nm and the right 
figure for A = 633nm. ^ le refractive index is m - 1.559 + i ■ 10“^ and m - 1.542 + i ■ 10“^ for A - 442nm and A - 633nm, 
respectively dOra 


in the paper, we use volume equivalent radii to represent the 
size of the particles, although the measured size distribution is 
obtained for projected-surface-area equivalent spheres. 

Fig .|3shows calculated polarization curves for distributions 
of hollow spheres with /max = 0.3,0.4 and 0.5 together with 
the laboratory measurements. The results for /max = 0.4 agree 
best with the measurements. This shape distribution reproduces 
the measured shape of the degree of linear polarization as a 
function of scattering angle reasonably well. The calculations 
show a feature at ~ 170° which may be caused by the fact that 
the particles in the hollow sphere distribution are still spherical 
in shape. 

In most applications the real size distribution of the dust 
grains is not known. In these cases one would like to extract 
the size distribution from the measured polarization curve. To 
do this we try to minimize the difference between the measure¬ 
ments and the calculations by fitting the size distribution. To 
find the least squares solution we have to minimize with re¬ 
spect to the size distribution. The;^f^ for this problem is defined 
by 




2 


Z 


LiW{ri) {Fu(9j,n)) 


2 


- PobsiOj) 


(7) 


In this equation {Fnk(0j, n)) is the n, kth element of the scatter¬ 
ing matrix of an ensemble of randomly oriented particles with 
a certain shape distribution and a volume equivalent radius r, 
at scattering angle 6j, w(r,) is the value of the relative num¬ 
ber size distribution at radius r, (these are the fitting parame¬ 
ters), PobsiOj) is the observed degree of linear polarization at 
scattering angle 6j and o-j is the error in PobsiOj)- The r, repre¬ 
sent different values in the size distribution, and the Oj represent 
the scattering angles at which measurements have been done. 
Unfortunately, Eq. 0 is not linear in the w{n), which makes 
minimization of very difficult, especially in the case when 


we want to sample the size distribution with high resolution, 
which means we would have a large number of fitting parame¬ 
ters. However, we can linearize the equation by minimizing 

-1:,' win) {{PniOj, n)) - PobsiOj) {PniOj, n )))" 

- . (8 

yj o-j 

where the yj are scaling factors which can be different for each 
6j. The right hand side of Eq. reduces to that of Eq. 0 when 

yj = ^ win) {puiOj,n)) ■ (9) 

i 

We can now iteratively solve Eq. 0 and 0 in the following 
way. 

1. As a first estimate take yj - 1 for all / 

2. Solve Eq 0 for the w(r,). 

3. The resulting values of win) are used to generate new yj 
from Eq. 0. 

4. Repeat steps 2 and 3 until convergence is reached, i.e. the 
yj do not change during one iteration. 

In all cases that we considered convergence is reached within 
15 iterations. This method is much faster than a nonlinear mini¬ 
mization algorithm and has the advantage that it is less sensitive 
to local minima. 

If one would use the above scheme without further consid¬ 
erations, the solution for the w(r,) that is found will in general 
not be a smooth function of the r, . Therefore, we add a regular¬ 
ization term to Eq. 0 which is defined as 

Tfi = ^ \win) - win-n)]^ , (10) 

i 

and minimize x’'^ + ^ Tr using the iterative scheme described 
above. The value of b is chosen in such a way that the con¬ 
tributions from the two t erms are balanced. Eor details on this 
regularization we refer to lPress et al.lill 9921) . 
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Fig. 9. The degree of linear polarization (upper panels) and the phase function Fi i (lower panels) of small quartz particles versus 
scattering angle, 6. All phase functions are normalized so that their values at 0 = 30° equals unity. Shown are the measurements 
(points) and the best fit to the polarization measurements using the distribution of hollow spheres with /max = 0.4 (solid curve). 
The left figures are for a wavelength A = 422 nm and the right figures for A = 633 nm. The refractive indices are as in Fig.|3 The 
curves for the two different wavelengths have been fitted independently and the resulting size distributions are plotted in Fig. cni 
The phase functions have not been fitted but are calculated using the size distribution derived from the fit to the polarization 
measurements. 


We have to avoid finding solutions with w(r,) < 0 or the 
solution with w(r,) = 0 for all i. To do this we use a linear 
least squares fitting procedure with linear equality and inequal¬ 
ity constraints. The constraints we use are 2, w(r,) = 1 and 
w(r;) > 0 for all i. 


From the derived w(r,) we can construct the relative 
projected-surface-area distribution S (log r). This size distribu¬ 
tion is defined in such a way that S (log r) d log r gives the rel¬ 
ative contribution to the total projected surface area of grains 
in the size range l og r to log r -i- d log r (see the appendix of 
l^lten et al.[Eoo3l . We have to note here that we use volume 
equivalent radii of the particles while the measured size distri¬ 
bution is based on a projected-surface-area equivalent radius. 
This might introduce a small difference in the size distributions 
when the particles deviate strongly from homogeneous spheres. 
The differences here will be minimal. 


The best fit results are obtained when we take /max = 0.4. 
The resulting polarization curves are shown in the upper pan¬ 


els of Fig. 13 Plotted in the lower panels of this figure are the 
phase functions Fn for the two different wavelengths. We see 
that the polarization curves are reproduced very well for both 
wavelengths. We have to note that for the other elements of 
the scattering matrix the agreement with the laboratory mea¬ 
surements is less convincing. For example depolarization and 
circular polarization effects cannot be explained using hollow 
spheres, because of the spherical symmetry of the particles. 
Also, as can be seen from Fig.|3 the calculated phase functions 
display an increase at large scattering angles which is absent 
in the measurements. The size distribution obtained by fitting 
the degree of linear polarization is plotted in Fig. cni together 
with the measured size distribution using laser diffraction. We 
see that the grain sizes contributing most to the total projected- 
surface-area of all three size distributions match very well. In 
general, we find that an extra contribution of very small grains 
is required with respect to the measured size distribution in 
order to reproduce the measured polarization curve. The laser 
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r (fim) 


Fig. 10. The size distribution, S (log r), of small quartz particles 
as measured by laser diffraction (solid curve) and as btted from 
the linear polarization curve at A - 442 nm (dashed curve) and 
at 4 = 663 nm (dotted curve). 


diffraction method used to measure the size distribution of the 
quartz particles in the laboratory is inaccurate for very small 
grains, which could be the reason why this small grain com¬ 
ponent is not found when this method is used. Another reason 
could be that the small grain component is used by the fitting 
procedure to simulate scattering by small scale structures in the 
particles itself. Considering the limited information used in the 
fitting procedure and the relatively simple grain shape distribu¬ 
tion used, the similarity is remarkable. This is an encouraging 
result for the use of the distribution of hollow spheres in inverse 
problems to obtain the size distribution from the degree of lin¬ 
ear polarization. We have to stress here that having a reasonable 
estimate of the value of the refractive index is important for ob¬ 
taining reliable results. 


5. Discussion 

In the previous sections it has been shown that the distribution 
of hollow spheres is a useful tool for calculating the optical 
properties of irregularly shaped particles. The question remains 
why this is the case. It is clear that the particles in astrophys- 
ical environments or in the laboratory are not hollow spheres. 
Therefore, there must be a fundamental characteristic of the 
particles that determines their optical properties. This charac¬ 
teristic on the one hand has to be important for the scattering 
and absorption properties and on the other hand it has to be 
clearly different for homogeneous spheres and hollow or non- 
spherical particles. One of the main things that discriminates a 
homogeneous sphere from other types of particles is its perfect 
symmetry. This perfect symmetry makes the interaction of an 
electromagnetic wave with this particle very sensitive to inter¬ 
ference effects. This will determine to a large extent the optical 
properties. Only small deviations from a homogeneous sphere 
are required to disturb these interference effects. This suggests 
that the first disturbance in the perfect homogeneous sphere is 
very important and is not very shape-dependent. 


6. Conclusions 


In this paper we examine the applicability of the statistical ap¬ 
proach for computations of optical properties of ensembles of 
irregularly shaped particles in random orientation when a dis¬ 
tribution of hollow spherical particles is used. 

The range of particle sizes for which we did numerical 
computations covers a critical region if we wish to study the 
combined size and shape effects on the absorption, extinction 
and scattering cross sections as functions of the wavelength. 
For smaller particles we are in the Rayleigh regime, where the 
shapes of the spectra are independent of particle size. For par¬ 
ticles larger than in the critical regime the effects of particle 
shape become negligible. We may thus conclude that, for all 
particle sizes, the emission and extinction spectra in the in¬ 
frared computed by using the distribution of hollow spheres 
display the same trends and overall spectral shape as a uniform 
distribution of spheroids. These spectra are significantly differ¬ 
ent from those obtained using homogeneous spheres. This sug¬ 
gests that the shape of these spectra is, to a certain degree, not 
sensitive to the exact particle shape as long as the deviations 
from homogeneous spherical particles are l arge enough. This 
result is in agreement with the conclusion of iMin et'af 112003 al) 
for very small particles and is extended here to particles of ar¬ 
bitrary size. 

The exact wavelength positions of observed spectral fea¬ 
tures are frequently used to identify the chemical composition 
of solid state particles. These wavelength positions also depend 
on particle size and shape. A reliable identification of solid state 
components can therefore only be made if the effects of both 
particle size and shape are taken into account. Calculations us¬ 
ing the distribution of hollow spheres allow such a detailed 
analysis. 

The degree of linear polarization for incident unpolarized 
light of hollow spheres with various volume fractions of the 
central vacuum inclusion was studied and compared to that of 
laboratory measurements at wavelengths of 442 and 633 nm. 
Good agreement was found with laboratory measurements for 
small, irregular quartz particles for scattering angles between 
5° and 173°. It was also shown that the particle size distribution 
can be derived by fitting the calculated degree of polarization 
to the laboratory data. With this procedure we found an over¬ 
abundance of small grains with respect to the measured size 
distribution. This can be due to the fact that the laser diffraction 
method used to measure the size distribution is inaccurate for 
small grains, or that the fitting procedure tries to simulate small 
scale structures that are present in the ’real’ particles but ab¬ 
sent in the hollow spheres. If the cause of the overabundance of 
small grains is the inaccuracy of laser diffraction when the par¬ 
ticles are small, fitting the degree of linear polarization might 
be a good alternative for obtaining information on the particle 
size distribution. 

Although our results imply that the optical properties of 
dust grains only contain limited information on the exact parti¬ 
cle shape, by using the distribution of hollow spheres we have 
a more reliable diagnostic to derive the mineralogy, size and 
total mass of these dust grains than when only homogeneous 
grains and CDE are used. We conclude that, because of its sim- 








12 


M. Min et al.: Modeling optical properties of cosmic dust grains using a distribution of hollow spheres 


plicity, low computational demand and success in reproducing 
the optical properties of irregular particles, the distribution of 
hollow spheres can be a powerful tool in studies of light scat¬ 
tering, absorption and emission by particles, and in particular 
in cases where large numbers of particle parameters need to be 
considered. 

Acknowledgements. We are grateful to L. B. F. M. Waters and 
H. Volten for enlightening conversations. We would like to thank 
M. I. Mishchenko and O. Munoz for valuable comments on an ear¬ 
lier version of this paper. 

References 

Aden, A. L. & Kerker, M. 1951, J. Appl. Phys., 22, 1242 
Bohren, C. F. & Huffman, D. R. 1983, Absorption and 
Scattering of Light by Small Particles (New York: Wiley) 
Bouwman, J., de Koter, A., Dominik, C., & Waters, L. B. F. M. 

2003, Astron. Astrophys., 401, 577 
Bouwman, J., Meeus, G., de Koter, A., et al. 2001, Astron. 
Astrophys., 375, 950 

Crovisier, J., Leech, K., Bockelee-Morvan, D., et al. 1997, 
Science, 275, 1904 

Dorschner, J., Begemann, B., Henning, T., Jaeger, C., & 
Mutschke, H. 1995, A&A, 300, 503 
Draine, B. T. 1988, ApJ, 333, 848 
Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 
Fabian, D., Henning, T., Jager, C., et al. 2001, Astron. 
Astrophys., 378, 228 

Fabian, D., Jager, C., Henning, T., Dorschner, J., & Mutschke, 
H. 2000, A&A, 364, 282 

Gray, D. E. 1963, American Institute of Physics (AIP). 
Handbook (New York: McGraw-Hill, 1963, 2nd ed., edited 
by Gray, Dwight E.) 

Greenberg, J. M. & Li, A. 1996, A&A, 309, 258 
Hanner, M. S., Lynch, D. K., & Russell, R. W. 1994, ApJ, 425, 
274 

Henning, T. & Stognienko, R. 1996, A&A, 311, 291 
Hildebrand, R. H. & Dragovan, M. 1995, ApJ, 450, 663 
Honda, M., Kataza, H., Okamoto, Y. K., et al. 2004, ApJ, 610, 
L49 

Hoogzaad, S. N., Molster, E. J., Dominik, C., et al. 2002, A&A, 
389, 547 

Hovenier, J. W. & van der Mee, C. V. M. 2000, 
in Light Scattering by Nonspherical Particles: Theory, 
Measurements, and Applications, ed. M. 1. Mishchenko, 
J. W. Hovenier, & L. D. Travis (San Diego: Academic Press), 
61-85 

Hovenier, J. W., Volten, H., Munoz, O., van der Zande, W. J., 
& Waters, L. B. 2003, Journal of Quantitative Spectroscopy 
and Radiative Transfer, 79, 741 

Jager, C., Dorschner, J., Mutschke, H., Posch, T., & Henning, 
T. 2003, A&A, 408, 193 
Jones, A. P. 1988, MNRAS, 234, 209 

Kahnert, E. M., Stamnes, J. J., & Stamnes, K. 2002, Journal of 
Quantitative Spectroscopy and Radiative Transfer, 74, 167 
Kahnert, M. 2003, Journal of Quantitative Spectroscopy and 
Radiative Transfer, 79, 775 


—. 2004, Journal of Quantitative Spectroscopy and Radiative 
Transfer, 85, 231 

Kemper, F, Vriend, W. J., & Tielens, A. G. G. M. 2004, ApJ, 
609, 826 

Kolokolova, L., Jockers, K., Gustafson, B. S., & Lichtenberg, 
G. 2001, J. Geophys. Res., 106, 10113 
Levasseur-Regourd, A.-C. 1999, Space Science Reviews, 90, 
163 

Li, A. & Draine, B. T. 2001, ApJ, 550, L213 

Li, A., Greenberg, J. M., & Zhao, G. 2002, MNRAS, 334, 840 

Mie, G. 1908, Ann Phys., 25, 377 

Min, M., Hovenier, J. W., & de Koter, A. 2003a, A&A, 404, 35 
—. 2003b, J. Quant. Spectrosc. Radiat. Transfer, 79-80, 939 
Mishchenko, M. L, Hovenier, J. W., & Travis, L. D., eds. 
2000, Light Scattering by Nonspherical Particles, Theory, 
Measurements and Applications (San Diego: Academic 
Press) 

Mishchenko, M. 1. & Travis, L. D. 1994, Applied Optics, 33, 
7206 

—. 1998, J. Quant. Spectrosc. Radiat. Transfer, 60, 309 
Mishchenko, M. L, Travis, L. D., & Lacis, A. A. 2002, 
Scattering, Absorption and Emission of Light by Small 
Particles (Cambridge: Cambridge University Press) 
Mishchenko, M. L, Travis, L. D., & Macke, A. 1996a, Applied 
Optics, 35, 4927 

Mishchenko, M. L, Travis, L. D., & Mackowski, D. W. 1996b, 
J. Quant. Spectrosc. Radiat. Transfer, 55, 535 
Molster, P. J., Waters, L. B. E. M., Trams, N. R., et al. 1999, 
A&A, 350, 163 

Munoz, O., Volten, H., de Haan, J. E, Vassen, W., & Hovenier, 
J. W. 2000, A&A, 360, 777 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Elannery, 
B. P. 1992, Numerical Recipes in fortran, 2nd edn. 
(Cambridge: Cambridge University Press) 

Purcell, E. M. 1969, ApJ, 158, 433 
Purcell, E. M. & Pennypacker, C. R. 1973, ApJ, 186, 705 
Servoin, J. L. & Piriou, B. 1973, Phys. Stat. Sol. (b), 55, 677 
Snow, T. P. & Witt, A. N. 1995, Science, 270, 1455 
Toon, O. B. & Ackerman, T. P. 1981, Appl. Opt., 20, 3657 
van Boekel, R., Waters, L. B. F M., Dominik, C., et al. 2003, 
A&A, 400, L21 

van de Hulst, H. C. 1957, Light Scattering by Small Particles 
(New York: Wiley) 

Volten, H., Munoz, O., Hovenier, J. W., et al. 2004, J. Quant. 

Spectrosc. Radiat. Transfer, in press 
Volten, H., Munoz, O., Rol, E., et al. 2001, J. Geophys. Res., 
106, 17,375 

Warren, J. L., Barret, R. A., Dodson, A. L., Watts, L. A., & 
Zolensky, M. E., eds. 1994, Cosmic Dust Catalog, Vol. 14, 
NASA Johnson Space Center, Houston 
Wriedt, T. 2002, Part. Part. Syst. Charact., 19, 256 



